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Abstract: We propose a new global and fully inclusive variable s^ in for determining the 
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mass scale of new particles in events with missing energy at hadron colliders. We define 
as the minimum center-of-mass parton level energy consistent with the measured values of the 
total calorimeter energy E and the total visible momentum P. We prove that for an arbitrary 
event, s l J^ n is simply given by the formula s l J^ n = \J E 2 - P% + y ^ E\ + Mf nv , where M inv is 
the total mass of all invisible particles produced in the event. We use tt production and several 
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supersymmetry examples to argue that the peak in the s^in distribution is correlated with 
the mass threshold of the parent particles originally produced in the event. This conjecture 
allows an estimate of the heavy superpartner mass scale (as a function of the LSP mass) in a 
completely general and model-independent way and without the need for any exclusive event 
reconstruction. In our SUSY examples of several multijet plus missing energy signals, the 
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accuracy of the mass measurement based on s„^ n is typically at the percent level, and never 
worse than 10%. After including the effects of initial state radiation and multiple parton 
interactions, the precision gets worse, but for heavy SUSY mass spectra remains ~ 10%. 
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i. 


Introduction 





The ongoing Run II of the Fermilab Tevatron and the imminent run of the Large Hadron 
Collider (LHC) at CERN are on the hunt for new physics beyond the Standard Model (BSM) 
at the TeV scale. Arguably the most compelling phenomenological evidence for BSM particles 
and interactions at the TeV scale is provided by the dark matter problem [1], whose solution 
requires new particles and interactions BSM. A typical particle dark matter candidate does 
not interact in the detector and can only manifest itself as missing energy. At hadron colliders, 
where the total center of mass energy in each event is unknown, the missing energy is inferred 
from the imbalance of the total transverse momentum of the detected visible particles, and 
is commonly referred to as "missing transverse energy" (MET). The dark matter problem 
therefore greatly motivates the study of MET signatures at the Tevatron and the LHC [2] . 

While the MET class of BSM signatures is probably the best motivated one from a 
theoretical point of view, it is also among the most challenging from an experimental point 
of view. On the one hand, to get a good MET measurement, one needs to have all detector 
components working properly, since the mismeasurement of any one single type of objects 
would introduce fake MET. In addition, there are complications from cosmics, pile-up, beam 
halo, noise, etc. Therefore, establishing a MET signal due to some new physics is a highly 
non-trivial task [2,3]. 

At the same time, interpreting a missing energy signal of new physics is quite challenging 
as well. The main stumbling block is the fact that we are missing some of the kinemati- 
cal information from each event, namely the energies and momenta of the missing invisible 
particles. What is worse, a priori we cannot be certain about the exact number of missing 
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Figure 1: The generic event topology under consideration in this paper. Black (red) lines correspond 
to SM (BSM) particles. The solid lines denote SM particles Xi, i = 1,2,... , n V i s , which are visible 
in the detector, e.g. jets, electrons, muons and photons. The SM particles may originate either from 
initial state radiation (ISR), or from the hard scattering and subsequent cascade decays (indicated 
with the green-shaded ellipse). The dashed lines denote neutral stable particles %i, i = 1, 2, . . . , rii nv , 
which are invisible in the detector. In general, the set of invisible particles consists of some number 
n x of BSM particles (indicated with the red dashed lines), as well as some number of 
SM neutrinos (denoted with the black dashed lines). The identities and the masses of the BSM 
invisible particles Xii (* = 1,2,..., n x ) do not necessarily have to be all the same, i.e. we allow for 
the simultaneous production of several different species of dark matter particles. The global event 
variables describing the visible particles are: the total energy E, the transverse components P x and 
P v and the longitudinal component P z of the total visible momentum P. The only experimentally 
available information regarding the invisible particles is the missing transverse momentum 



particles in the event, or their identity, e.g. are they SM neutrinos, new BSM dark matter 
particles, or some combination of both? These difficulties are illustrated in Fig. |], where we 
show the generic topology of the missing energy events that we are considering in this paper. 
As can be seen from the figure, we are imagining a completely general setup - each event 
will contain a certain number n v i s of Standard Model (SM) particles Xi, i = 1,2,... ,n v i s , 
which are visible in the detector, i.e. their energies and momenta are in principle measured. 
Examples of such visible SM particles are the basic reconstructed objects, e.g. jets, photons, 
electrons and muons. The visible particles Xi are denoted in Fig. [l] with solid black lines 
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and may originate either from initial state radiation (ISR), or from the hard scattering and 
subsequent cascade decays (indicated with the green-shaded ellipse). On the other hand, 
the missing energy fix (or more appropriately, the missing transverse momentum fx) will 
arise from a certain number rii nv of stable neutral particles %j, i = 1, 2, . . . , rii nv , which are 
invisible in the detector. In general, the set of invisible particles in any event will consist 
of a certain number n x of BSM particles (indicated with the red dashed lines), as well as a 
certain number n u = rii nv — n x of SM neutrinos (denoted with the black dashed lines). The 
missing energy measurement alone does not tell us the number rii nv of missing particles, nor 
how many of them are neutrinos and how many are BSM (dark matter) particles. Notice 
that in this general setup the identities and the masses rrii of the BSM invisible particles Xii 
(i = 1,2,..., n x ) do not necessarily have to be all the same, i.e. we allow for the simultaneous 
production of several different species of dark matter particles [4-7]. On the other hand, we 
shall always take the neutrino masses to be zero 

mi = 0, for i = n x + l,n x + 2, ...,m nv . (1.1) 

Most previous studies of MET signatures have assumed a particular BSM scenario and 
investigated its consequences in a rather model-dependent setup. The results from those 
studies would seem to indicate that in order to make any progress towards determining what 
kind of new physics is being discovered, and in particular towards mass and spin measure- 
ments, one must attempt at least some partial reconstruction of the events, by assuming a 
particular production mechanism, and then identifying the decay products from a suitable 
decay chain [8-56]. In doing so, one inevitably encounters a combinatorial problem whose 
severity depends on the new physics model and the type of discovery signature. For example, 
complex event topologies with a large number n v i s of visible particles, and/or a large number 
of jets but few or no leptons, will be rather difficult to decipher, especially in the early data. 
Therefore, it is fair to ask whether one can say something about the newly discovered physics 
and in particular about its mass scale, using only inclusive and global 1 event variables, before 
attempting any event reconstruction. 

In this paper, therefore, we shall concentrate on the most general topology exhibited 
in Fig. |] and we shall make no further assumptions about the underlying event structure. 
For example, we shall not specify anything about the production mechanism. In particular, 
we shall not make the usual assumption that the BSM particles are pair produced and, 
consequently, that there are two and only two BSM decay chains resulting in n x = 2 identical 
dark matter particles with equal masses m\ = ni2- Accordingly, we shall not make any 
attempt to group the observed SM objects JQ, i = 1,2,..., n v i s , into subsets corresponding 
to individual decay chains. Furthermore, we shall in principle allow for the presence of SM 
neutrinos which could contribute towards the measured MET. In this sense our approach will 
be completely general and model- independent. 

Given this very general setup, our first goal will be to define a global event variable 
which is sensitive to the mass scale of the particles that were originally produced in the 

1 Here and throughout the paper, we use the term "global" from an experimentalist's point of view. Strictly 
speaking, the detectors are not fully hermetic, hence no variable can be truly global in the theorist's sense. 
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event of Fig. [T], or more generally, to the typical energy scale of the event. Since we are 
not attempting any event reconstruction, this variable should be defined only in terms of 
the global event variables describing the visible particles X{, namely, the total energy E 
in the event, the transverse components P x and P y and the longitudinal component P z of 
the total visible momentum P in the event. In the same spirit, the only experimentally 
available information regarding the invisible particles that we are allowed to use is the missing 
transverse momentum Pt (see Fig. ||). Of course, the missing transverse momentum Pt is 
related to the transverse components P x and P y of the total visible momentum P as 

f T = - (P x e x + PySy) = -P T , (1.2) 

so that we can use Pt and Pt = P x &x + Py&y interchangingly. Then, the commonly used 
missing energy Et is nothing but the magnitude Pt of the measured missing momentum Pt- 



E T ^P T = P T = ^P2 + P2 . (1.3) 
The main idea of this paper is to propose a new global and. inclusive variable Smm 

defined 

as follows. s m i n is simply the minimum value of the parton- level Mandelstam variable s which 
is consistent with the observed set of E, P z and Pt in a given event 2 . Correspondingly, its 
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square root s', is the minimum parton level center-of-mass energy, which is required in 
order to explain the observed values of E, P z and Et- Our main result, derived below in 
Section is the relation expressing the so defined s l J^ n in terms of the measured global and 
inclusive quantities E, P z and Et- In Section [2] we shall prove that sl^f is always given by 
the formula 

a 1 / 2 c\/r. \ — . Fv2 D2 1 . / 77-2 1 71^-2 



t n (M inv ) = VW^P! + sjE% + Ml v , (1.4) 

where the mass parameter Mi nv is nothing but the total mass of all invisible particles in the 
event: 



riini, n 



X 



M inv = nii = VJ rm , (1.5) 

8=1 1=1 

and the second equality follows from the assumption of vanishing neutrino masses (|1 . 1|) . 

As can be seen from its defining equation (0), the variable s^ 2 n is actually a function of 
the unknown mass parameter Mi nv . This is the price that we will have to pay for the model- 
independence of our setup. This situation is very similar to the case of the Cambridge Mt2 
variable [9,14,34-37,45,52,53] and its various cousins [33,38,39,41,42,46,48,50,54,55], which 
are also defined in terms of the unknown test mass of a missing BSM particle. However, 
the Cambridge Mt2 variable is a much more model-dependent quantity, since it requires 
the identification of two separate decay chains in the events. Furthermore, in some special 
cases (more precisely, those of M^' n ' n ^ in the language of [54]) Mt2 is essentially a purely 
transverse quantity, and in this sense would not make full use of all of the available information 



2 In what follows, instead of J?t we choose to use the more ubiquitous $t, since the two are essentially the 



same, see (1.3) 
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in the event. In contrast, our variable s', is denned in a fully inclusive manner, and uses 
the longitudinal event information as well. 

After deriving our main result (1.4) in Section B we devote the rest of the paper to 
studies of its properties. For example, in Section [3| we shall compare s l J^ n to some other 
global and inclusive variables which have been considered as measures of the mass scale of 
the new particles: Ht [12], the total visible invariant mass M [2], the missing transverse 
energy Et, the total energy E, and the total transverse energy Et in the event. We shall use 
several examples from SM ti production, as well as supersymmetry (SUSY), to demonstrate 
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that among all those possibilities, the variable s„^ n is the one which is best correlated with 
the mass scale of the produced particles, even when we conservatively set the unknown mass 
parameter Mi nv to zero. In Section ||] we shall investigate the dependence of the s l J^ n variable 
on the a priori unknown mass parameter Mi nv , using conventional SUSY pair-production for 
illustration. We shall find a very interesting result: when the parameter M{ nv happens to be 

1/2 

equal to its true value, the peak in the distribution is surprisingly close to the SUSY 
mass threshold. This correlation persists even when the two SUSY particles produced in the 
hard scattering are very different, for example, in associated gluino-LSP production. This 
observation opens up the possibility of a new, all inclusive and completely model-independent 
measurement of the mass scale of the new (parent) particles produced in the event: we 

,1/2 

simply read off the location of the peak in the distribution, and interpret it as the 
mass threshold of the parent particles. Because of the intrinsic dependence on the unknown 
mass parameter Mi nv , the method only provides a relation between the mass of the parent 
particle and the mass of the dark matter particle, just like the method of the Cambridge Mt2 
variable [9]. However, unlike the Mt2 endpoint measurements, our measurement is based on 
an all-inclusive global variable, and does not require any event reconstruction at all. It is 
worth noting that since we are correlating a physics parameter to the peak, rather than the 
endpoint of an observed distribution, our measurement will be less prone to errors due to finite 
statistics, detector resolution, finite width effects etc., which represents another important 

1/2 

advantage of the s^ in variable. The accuracy of our new mass measurement method is 
investigated quantitatively in Sections [p and ^. Our discussion in Sections ||, || and [|, while 
demonstrating the usefullness of the sj? variable, will be limited to an ideal case, where 
the effects from initial state radiation (ISR), multiple parton interactions (MPI) and pile-up 
are negligible. In Section || we investigate the adverse effect of those latter factors on the 
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Smin measurement in a realistic experimental environment and discuss different approaches 
for minimizing their impact. In Section ^ we summarize our main points and conclude. 



2. Derivation of s^ in 

In this section we shall derive the general formula ( |1.4| ) advertised in the Introduction. Before 
we begin, let us introduce some notation. We shall denote the three-momenta of the invisible 
particles Xii * = 1> 2, . . . , rii nv , with pi, or in components Pi x , pi y and pi z . As usual, we 
choose the z-axis along the beam direction, so that pi x and pi y are the components of the 
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transverse momentum pix- As already mentioned in the Introduction, the masses of the 
invisible particles will be denoted by mi. 

Our starting point will be the expression for the parton-level Mandelstam variable s for 
the event depicted in Fig. [I]: 



2 



i=l / V i=l 



The invisible particle momenta pi are not measured and are therefore unknown. However, 
they are subject to the missing energy constraint: 

Y^PiT = Ft = -Pt , (2-2) 
i=i 

which causes the second term in fl2.1|) to vanish and we arrive at a simpler version of fl2.1| ) 

E + Y.^ m2 i+P?T+Pl \ -(ft + Eft.J • ( 2 - 3 ) 

We see that the expression for s is a function of a total of 3rii nv variables pi which are subject 
to the 2 constraints ( |2.2j ). Given that we are missing so much information about the missing 
momenta pi, it is clear that there is no hope of determining s exactly from experiment, and the 
best one can do is to use some kind of an approximation for it. For example, Ref. [52] recently 
proposed to approximate the real values of the missing momenta pi with the values that 
determine the event Mt2 variable. However, constructing any Mt2 variable requires one to 
make certain model-dependent assumptions about the underlying topology of the event, and 
furthermore, for very complex events, with large n v i s , the associated combinatorial problem 
will become quite severe. Therefore, here we shall use a different, more model-independent 
approach. The key is to realize that the function s has an absolute global minimum s m i n ^ when 
considered as a function of the unknown variables pi. Therefore, we choose to approximate 
the real values of the missing momenta with the values corresponding to the global minimum 



s m i n . The minimization of the function (2.3) with respect to the variables pi, subject to the 



constraint fl2.2| ), is rather straightforward. The global minimum is obtained for 

PiT = T7^ Ft, (2.4) 



where the parameter 




(2.5) 



(2.6) 



- 6 - 



was already defined in ( |1.5[ ) and represents the total mass of all invisible particles in the event. 
Since the neutrinos are massless, Mi nv only counts the masses of the BSM invisible particles 
which are present in the event. Substituting ( |2.4| ) and into ( |2.3| ) and simplifying, we 
get the minimum value s m i n of the function ( |2.3| ) to be 

Smin{M inv ) = LfE* ~ P* + J?* + MjJ . (2.7) 

Since the right-hand side is a complete square, it is convenient to take the square root of 
both sides and consider instead 

^L(M mv ) = y/W=Pl+yJp*+Ml v , (2.8) 

which can be equivalently rewritten in terms of the missing energy Et as 

€UMinv) = y r E^Pl + y/p%+M? nv , (2.9) 

completing the proof of ( |1.4j) . 

A few comments regarding the variable s 1 J^ n defined in (^1]) are in order. Perhaps the 
most striking feature of s l J^ n is its simplicity: the result ( |2.9| ) holds for completely general 

1 /2 

types of events, with any number and/or types of missing particles. Clearly, s^ in itself is both 
a global and an inclusive variable, since it is defined in terms of the global and inclusive event 
quantities E, P z and Et, which do not require any explicit event reconstruction. It is easy to 
see that the expression ( p.9| ) is invariant under longitudinal boosts, since it depends on the 
quantities E 2 — P 2 2 , Et and Mi nv , all three of which are invariant under such boosts. Also 

1/2 

notice that s^ in has units of energy and thus provides some measure of the energy scale in 
the event, and can be directly compared to other popular energy-scale variables (see Section |3] 
below). In the remainder of this paper we shall investigate in more detail the properties of 
the new variable (|2.S|). 

~l/2 

3. Comparison between s^ in and other global inclusive variables 

The immediate question after the discovery of a MET signal of new physics at the Tevatron 
or LHC, will be: "What is the energy scale of the new physics?". We shall now argue that 
our global inclusive variable s 1 J^ n from ( |2.9D provides a first, relatively quick answer to this 
question, which will turn out to be surprisingly accurate, given that we are not attempting 
any event reconstruction or modelling of the new physics. Of course, one might do better 
by considering exclusive signatures and applying the usual tricks for mass measurements, 
but chances are that this will require some time. It is therefore worth investigating how 
much information one can get from totally inclusive measurements like (|2,9| ) which should be 
available from very early on. 

To set up the subsequent discussion, let us introduce the different global variables from 
Fig. H which will be experimentally accessible. The total visible energy E is simply 

E = J2^ a , (3.1) 

a 
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where we use the index a to label the calorimeter towers, and E a is the energy deposit in 
the a tower 3 . As usual, since muons do not deposit significantly in the calorimeters, the 
measured E a should first be corrected for the energy of any muons which might be present 
in the event and happen to pass through the corresponding tower a. The three components 
of the total visible momentum P are 

P x = ^2 E a sin 9 a cos ip a , (3.2) 

a 

P y = Ea sin 9 a sin (p a , (3.3) 

a 

P z =J2 E aCOS9 a , (3.4) 

a 

where 9 a and ip a are correspondingly the azimuthal and polar angular coordinates of the a 
calorimeter tower. The total transverse energy Ey? is 

E T = J2 E aSin9 a , (3.5) 



while the missing transverse energy Et was already defined in (1.3). 

We are now in a position to introduce the variable Ht which is commonly used throughout 
the literature, yet, quite surprisingly, there is no universally accepted definition for it. The 
idea behind Ht is to add up the transverse energies of various objects in the event, including 
the missing energy ( |1.3| ). While the idea is rather straightforward, there are large variations 
when it comes to its implementation. For example, one issue is whether one should use only 
reconstructed objects or simply sum over all calorimeter towers as we have been doing here 
so far. The former method has the advantage that it would tend to reduce pollution from 
the underlying event, noise, etc. On the other hand, it would introduce dependence on the 
jet reconstruction algorithm, the ID cuts, etc. Those subtleties are avoided in the second 
method, which defines a purely calorimeter based Ht- There are other possible variations in 
the definition of Ht, for example, whether one includes all jets, or just the top 4 in pt [12], 
whether or not one includes the leptons in the sum, etc. For the purposes of this paper, we 
do not need to go into such details, and we shall simply use a calorimeter-based, all inclusive 
Ht definition as 

Ht = Et~\- Et • (3-6) 
Finally, we shall also consider the total visible mass in the event [2] 



M = y/E* - Pi - PI - Pi = X /E*- P T - PI . (3.7) 

n 

alternatively written in a more symmetric form as 



Note that in terms of the visible mass M just introduced, our s l J^ n variable fl2,9| ) can be 



S^niMinv) = \/E T + + J El + Mf nv . (3.8) 



3 We ignore the difference in the segmentation of the hadronic and electromagnetic calorimeters, and for 
E a simply add up the HCAL and ECAL energy deposits. 
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Figure 2: Unit-normalized distributions of the various energy scale variables Ei introduced in Sec- 
tion^: E (blue), Et (cyan), Et (magenta), Ht (green), M (red) and sHf n (0) (black); in (a) single- 
lepton and (b) dilepton ti events. The dotted (yellow-shaded) histograms are identical in panels (a) 
and (b) and show the true s 1 / 2 distribution. 

We are now ready to contrast the so denned global inclusive variables E, Et, Et, Ht and 
M to our variable §Hf n denned in (|2.9|), Since sl[f n (M inv ) depends on the a priori unknown 
invisible mass parameter M{ nv , first we need to decide what to do about the M{ nv dependence 
in (^1]). In the remainder of this section, we shall adopt a most conservative approach: we 
will simply set Mi nv = and consider the variable 

B^{0) = y/&=P*+ . (3.9) 

This choice is indeed very conservative: for SM processes, where the missing energy is due 
to neutrinos, this would be the proper variable to use anyway. On the other hand, for 
BSM processes with massive invisible particles, at this point we are lacking the necessary 
information to make a more informed choice. We shall postpone our quantitative discussion 
of the Mi nv dependence in ( |2.9| ) until the next section |j. 

We shall illustrate our comparisons with specific examples, illustrated in Figs, g ^ and 
^|. In each case, we shall plot the six different global inclusive variables Ei introduced so far, 
with the following color scheme: in Figs. we shall plot the calorimeter energy E (|3.l| ) 
with blue lines, the missing transverse energy Et (|1.3| ) with cyan lines, the total transverse 
energy Et ( |3.5[) with magenta lines, the Ht variable (|3.6| ) with green lines, the total visible 
mass M (|3.7|) with red lines, and finally, our s^L(O) variable ( |3.9|) with solid black lines. All 
numerical results shown here have been obtained with PYTHIA 4 [57] and the PGS detector 
simulation package [58]. As our first example, shown in Fig. |^, we choose ti production at the 
LHC (the corresponding data from the Tevatron already exists, so the same comparison can 
also be made directly with CDF and DO data as well). In Fig. ^(a) (Fig. ^(b)) we show our 

4 For simplicity, for the numerical results shown in this and the next two sections, we turned off ISR and 

1/2 

MP1 in PYTHIA, which allows us to better illustrate and subsequently explain the salient features of s^ in . 
The ISR and MPI effects will be studied later in Section |(| 



- 9 - 



results for the semi-leptonic (dilepton) channel. The dilepton tt sample is rather similar to 
a hypothetical new physics signal due to dark matter particle production: each event has a 
certain amount of missing energy, which is due to two invisible particles escaping the detector. 

In each panel of Fig. ^ the dotted (yellow-shaded) histogram shows the true s 1 / 2 dis- 
tribution, which is the one we would ideally want to measure. However, due to the missing 
neutrinos, s 1//2 is not directly observable, unless we make some further assumptions and at- 
tempt some kinematical event reconstruction. Therefore we concentrate on the remaining 
distributions shown in Fig. ^, which are immediately and directly observable. In particular, 
we shall pose the question, which among the various distributions exhibited in Fig. [2] seems to 
be the best approximation to the true s 1//2 distribution. A quick glance at Fig. || reveals that 
the variable which comes closest to the true s 1 / 2 is precisely our variable sllf n {0) defined in 
( |3.9| ). As for the rest, we see that the missing transverse energy Et is a very poor estimator 
of the energy scale of the events, while Et, Ht and M are doing a little bit better, yet are 
still quite far off. As can be expected from its definition (|3.6|), Ht is always somewhat larger 
than Et, while Ht and M are rather similar, with Ht (M) doing better for the dilepton 
(semi-leptonic) case. Finally, the total energy E is relatively close to the true s 1 / 2 distribu- 
tion, but is quite broad in both Figs. |2|(a) and §(b). In contrast, the s~Hf n (0) distribution is 
quite sharp, and is thus a better indicator of the relevant energy scale. 

Let us now take a closer look at the two s 1 / 2 distributions in each panel of Fig. |2[ 

1 Si 

Since s^ in was defined through a minimization procedure, it is clear that it will always 
underestimate the true s 1 / 2 . Fig. || quantifies the amount of this underestimation for the 
case of tt events. We see that sllf n (0) is tracking the true s 1 / 2 quite well for the case of 
semi-leptonic tt events in Fig. |2|(a). This could have been expected on very general grounds: 
for semi-leptonic events, we are missing a single neutrino, whose transverse momentum is 
actually measured through ftr, so that the only mistake we are making in approximating 
s 1 / 2 sHf n (0) is due to the unknown longitudinal component p\ z . In the case of dilepton 
events, however, there are two missing neutrinos, and thus more unknown degrees of freedom 



which we have to fix rather ad hoc according to our prescription (2.4, 2.5). The resulting 
error is larger and leads to a larger displacement between the true s 1 / 2 distribution and its 
sllf n (0) approximation, as can be seen in Fig. §(b). 

In the case of ti illustrated in Fig. || the missing energy arises from massless SM neutrinos, 
so that the approximation Mj nt) = is well justified. Let us now consider a situation where 
the observed missing energy signal is due to massive neutral stable particles, as opposed 
to SM neutrinos. The prototypical example of this sort is low energy supersymmetry with 
conserved i?-parity, and this is what we shall use for our next two examples as well. Each 
SUSY event will be initiated by the pair-production of two superpartners, which will then 
cascade decay to the lightest super symmetric particle (LSP), which we shall assume to be 
the lightest neutralino x?- Since there are two SUSY cascades per event, there will be two 
LSP particles in the final state, so that 

n inv = n x = 2 . (3.10) 
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Figure 3: The same as Fig. |2| but for gluino pair production events with (a) 2-jet gluino decays as in 
(gjgj and (b) 4-jet gluino decays as in ( |jg ). The SUSY masses are hxed as follows: m^o = 100 GeV, 



m^o — 200 GeV and rrig = 600 GeV. In addition to the variables shown in Fig. g, here we also plot the 
*mm(^ m x) distribution (dotted line) with the correct value of the invisible mass Mi nv = 2m x = 2m x o . 



Furthermore, since the two LSPs are identical, we also have 

mi = ?7i2 = Tn x , (3-H) 
i.e. in what follows we shall denote the true LSP mass with m x . From (|1.5| ), ( |3.10 ) and 



(3.11) it follows that the true total invisible mass in any SUSY event is simply 

M inv = 2m x . (3.12) 

However, the true LSP mass m x is a priori unknown, therefore, when we construct our 
variable 

$LU. M ™) = s]il(2m x ) (3.13) 

for the SUSY examples, we will have to make a guess for the value of the LSP mass m x . We 
shall denote this trial value by fh x , in order to distinguish it from the true LSP mass m x . 
This situation is reminiscent of the case of the Cambridge Mt2 variable [9] , where in order to 
construct the Mt2 variable itself, one must first choose a test value for the LSP mass. Our 
notation here is consistent with the notation for M^2 used in [54]. 

We are now ready to describe our SUSY examples. For our study we will choose a 
rather difficult signature — jets plus for which all other proposed methods for mass 
determination are bound to face significant challenges. For concreteness, we consider gluino 
production, followed by a gluino decay to jets and a neutralino. In Fig. |3| we consider gluino 
pair-production (gg), while in Fig. |4| we show results for associated gluino-LSP production 
(gXi) ■ in addition, we consider two different possibilities for the gluino decays. The first 
case, shown in Figs. [|(a) and Figs. |](a), has the gluino decaying directly to the LSP: 

9 - jjXi , (3-14) 
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so that the gluino pair-production events in Fig. ||(a) have 4 jets and missing energy, while 
the associated gluino-LSP production events in Fig. |||(a) have two jets and missing energy. 
In the second case, presented in Figs. ||(b) and Figs. f|(b), we forced the gluino to always 
decay to X21 w hich in turn decays via a 3-body decay to 2 jets and the LSP: 

9 -> 33x1 ~* 3333X1 ■ (3-15) 

As a result, the gluino pair-production events in Fig. |3|(b) will exhibit 8 jets and missing 
energy, while the associated gluino-LSP production events in Fig. ||](b) will have four jets and 
missing energy. Of course, the actual number of reconstructed jets in such events may be even 
higher, due to the effects of initial state radiation (ISR) and/or jet fragmentation. In any 
case, such multijet events will be very challenging for any exclusive reconstruction method, 
therefore it is interesting to see what we can learn about them from the global inclusive 
variables discussed here. 

For concreteness, in what follows we shall always fix the relevant SUSY masses according 
to the approximate gaugino unification relation 

rrig = 3m^o = 6m^o = 6m x , (3.16) 



and since we assume three-body decays in ( 3.14 ) and ( |3.15 ), we do not need to specify the 



SUSY scalar mass parameters, which can be taken to be very large. In addition, as implied 
by ( 3.16| ), we imagine that the lightest two neutralinos are gaugino- like, so that we do not 



have to specify the higgsino mass parameter either, and it can be taken to be very large as 
well. 

Fig. H shows our results for the different global inclusive variables introduced earlier, for 
the case of gluino pair-production. All in all, the outcome is not too different from what 
we found previously in Fig. |2| for the ti case: when it comes to approximating the true s 1 / 2 

1/2 

distribution, the missing energy $t does the worst, our variable s^ in (0) does the best, and 
all other remaining variables are somewhere in between those two extremes. This time, in 
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Fig. |3] we also plot one "cheater" distribution, namely s m j„(2m x ), where we have used the 
correct value of the invisible mass Mi nv = 2m x = 2m^o . It demonstrates that knowing the 

1 /2 

actual value of the LSP mass helps (since ^ in (2m x ) gets closer to the truth), but is not 
crucial: the quantity sHf n (0) still does surprisingly well in approximating the true s 1 / 2 . 

Notice that when the missing energy in the data is due to massive BSM particles, there are 
two sources of error in approximating s 1 / 2 ~ Smi n (®)i eac h leading to an underestimation. By 
comparing the three different types of s 1 / 2 distributions shown in each panel of Fig. ||, one can 
see quantitatively the effect of each source. First, when we take the minimum possible value 
of s 1 / 2 in ( |2.3D , we are underestimating by a certain amount, which can be seen by comparing 
the "cheater" distribution s^ n (2m x ) (dotted line) to the s 1 / 2 truth (yellow shaded). Second, 
as we do not know a priori the LSP mass, we take conservatively M{ nv = 0, which leads to 

1 /2 

a further underestimation, as evidenced by the difference between the s„ in (0) distribution 

1 /2 

(solid line) and its "cheater" version •s r) (j„(2m x ). In spite of those two undesirable effects, the 

,1/2 

5^(0) approximation that we end up with is still surprisingly close to the real one, and is 
certainly the best approximation among the variables we are considering. 

The common thread in our first two examples shown in Figs. [2| and [| was that the events 
were symmetric, i.e. we produce the same type of particles, which then decay identically 
on each side of the event. As our last example, we shall consider an extreme version of an 
asymmetric event, namely one where all visible particles come from the same side of the event, 
i.e. from a single decay chain. The process of associated gluino-LSP production is exactly of 
this type - all jets arise from the decay chain of a single gluino, which is recoiling against an 
LSP. The topology of these events is very different from the events considered earlier in Figs. |2| 
and 0. Nevertheless, as seen in Fig. 01, we find very similar results. In particular, among all 
the different global inclusive variables that we are considering, the quantity s (0) is still 
the one closest to the true s 1//2 distribution. 



,1/2 

4. Dependence of s„( iri on the unknown masses of invisible particles 

In the previous Section |3] we demonstrated the advantage of s^ in in comparison to the other 
commonly used global inclusive event variables. From now on we shall therefore focus our 

1/2 

discussion entirely on s^ in and its properties. In this Section we shall investigate in more 

1/2 

detail the dependence of on the (a priori unknown) masses of the invisible particles which 
are causing the observed missing energy signal. Then in the next Section S we shall use these 

1/2 

results to correlate the observed s^ in distribution to the masses of the parent particles which 
were originally produced in the event. 

Recall that in the three examples from the previous section, we always conservatively 

1/2 

chose the invisible mass to be zero: Mi nv = and we correspondingly considered s^ in (0). 
This choice is actually a good starting point in studying any missing energy signature by 

1/2 

means of s^ in {Mi nv ). The assumption of Mi nv = is precisely what one would do if one 
were to assume that the missing energy is simply due to SM neutrinos, as opposed to some 
new physics. However, if the observed missing energy signal is in excess of the expected SM 



- 13 - 



> 2 

a 



z 

■a 



ii-4j+^ T (m x =100 GeV) (a) 
m = 100,200,300,400,500 




) 1500 i 

s^ 2 n (2m x ) (GeV) 



2500 



> 



Z 



0-3 
7 
5 



ii-4j+^ T (m x =200 GeV) (b) 

= 0,100,200,300,400,500 




2000 2500 3000 

si£(2m x ) (GeV) 



> 



>1 



z 

■d 



££-4j+£ T (m.^300 GeV) (c) 
i y = 0,200,300,400,500,600 




3000 



3500 4000 

s^(2m„) (GeV) 



> 



>1 



z 

T3 



ii-4j+^ T (m x =400 GeV) (d) 
m x = 200,400,600,800,1000 

1/2 




4500 



5000 



si£(2m x ) (GeV) 



6000 



Figure 5: Unit-normalized distributions of the s^^Minv) variable for several different SUSY mass 
spectra: (a) m^o = 100 GeV, (b) m^o = 200 GeV, (c) m^o = 300 GeV, and (d) m$> = 400 GeV. 
The remaining masses are fixed according to ( 3.ldQ . We consider gluino pair-production events with 



1/2 

2-jet gluino decays as in ( |3.l4 K In each panel, we plot the s^ n (M inv ) = s rj ( in (2TO x ) distributions 
for several representative values of the trial LSP mass m x as shown. The color scheme is such that 
the black histogram is always the case where we happen to use the correct value of the LSP mass: 
fh x = m x . The dotted (yellow-shaded) histogram gives the true s 1 / 2 distribution. 



backgrounds, then an alternative, BSM explanation for those events must be sought. In that 
case, we would not know the mass of the invisible particles, and we would have to make a 
guess. Our main goal in this section is to study numerically the effect of this guess. Our 
philosophy will be to revisit the SUSY examples from Section [3| and simply vary the test mass 
m x of the invisible particles (the LSPs). Since the two LSPs are identical (see eq. ( 3.11| )), we 
will take their test masses to be the same as well. 

Our results are presented in Figs. [5], || and [?| In Figs. |5| and || we consider gluino pair 
production. In Fig. |5| each gluino decays to 2 jets as in ( |3.14| ), while in Fig. ^ each gluino 
decays to 4 jets as in ( |3.15 ). Then in Fig. |7| we consider asymmetric events of associated 
gluino-LSP production, where the single gluino decays to 4 jets as in ( 3.15| ). In each figure, 
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Figure 6: Tie same as Fig. 0, but for 4-jet gluino decays as in (3.15) 



we consider four different study points, defined through the value of the true LSP mass m x . 
In all three Figs. H-0, panels (a) correspond to m x = 100 GeV, panels (b) have m x = 200 GeV, 
panels (c) have m x = 300 GeV, while in panels (d) m x = 400 GeV. As before, the remaining 
masses rrig and m^o are always fixed according to the approximate gaugino unification relation 
(3.16). Each panel in Figs. p|-f7| exhibits the true s 1 / 2 distribution (yellow-shaded histogram), 

~l/2 

and the corresponding s' in (2m x ) distributions for several representative values of the test 

,.1/2 

LSP mass m x . Each s '■ curve is both color coded and labelled by its corresponding value 
of rh x . Our color scheme is such that the s ' 4 histogram in black is the one where we happen 
to use the correct value of the LSP mass, i.e. when rh x = m x . 

The qualitative behavior seen in Figs. ^-0 is more or less as expected: the s]f^ n (2rh x ) 
distributions shift to higher energy scales, as we increase the value of the test mass m x . This 
can be easily understood from the definition ( |2 .9| ) of the s l J^ n {Mi nv ) variable: for any given 

1/2 

set of E, P z and Et values, s^^iM^) is a monotonically increasing function of M\ nv . The 
shifts observed in Figs. ||-|7] also make perfect physical sense: obviously, one needs more energy 
in order to produce heavier invisible particles. 
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Figure 7: The same as Fig. [g, but for events of associated gluino-LSP production (gxi)- 



Let us now concentrate on the quantitative aspects of Figs. H-0. Upon careful inspection 
of the three figures, we notice that when the test mass m x is equal to the true mass m x 

1/2 

(i.e. for the black colored histograms), the corresponding distribution s^ in (2m x ) peaks very 
close to the true s 1 / 2 threshold [s l ^ 2 \ hr - As usual, we define the threshold (s 1//2 ) tftr as the 
value where the true s 1 / 2 distribution (yellow shaded histogram) sharply turns on. This 
observation is potentially extremely important, since the threshold (•s 1//2 ) i?tr is simply related 
to the masses of the two particles which were originally produced in the event. For example, 
for the gluino pair production events in Figs. ^ and |6| the threshold is given by 



s 



1 / 2 \ = 2m- g = I2m x , (4.1) 



thr 



where the second equality is valid only under the gaugino unification assumption ( p,16| ). 
Similarly, in the case of associated gluino-LSP production in Fig. [?], the threshold is given by 

s 1 / 2 ) _ = m- g + m^o = 7m x , (4.2) 



/ thr 
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where once again the second equality is due to our assumption ( 3.1 6| ) . It is easy to verify that 



in all three figures ^, || and[^, the s 1//2 thresholds (i.e. the sharp turn-ons in the yellow-shaded 
distributions) always occur at the locations predicted in eqs. (4.1) and Q4.2p. 



Let us now introduce one last piece of notation. In what follows we shall use the notation 

(slSniMinv)) h (4.3) 
V / peak 

1/2 

to denote the particular value of s^ in where we find the peak of the distributions 

dN(s]il{M inv )) 
uo rmn 



= 0. (4.5) 

/ peak 



which are plotted in Figs. In other words, 

d dN{s^(M inv )) 



ds 1 ' 2 ds 1/2 

mm mm 



s ' ■ = 

mm 



With those conventions, we can now formulate our empirical observation above as 

(* 1/a ), ft - (4 7 l(2-x)) k ■ (4-6) 

V / thr \ / peak 

The last equation is one of the main results in this paper. While we were not able to derive it 
in a strict mathematical sense, it is nevertheless supported by our numerical results shown in 
Figs. [5||7[ We also checked many other SUSY examples, where we used different mass spectra 
and different production processes and decays. We found that in all cases the approximate 



relation (4.6) still holds. Fig. || quantifies this statement for the two previously considered 
processes of gluino pair production and associated gluino-LSP production, where the gluinos 
are forced to decay either to 2 jets as in ( fug) or to 4 jets as m fl3~T5l) . In the fi gure we 



compare the following three quantities, all of which are related in one way or another to the 
energy scale s 1//2 of the events: 

• (s 1 ^ 2 ) ave '■ this is the average of the true s 1//2 distribution (the one shown in the previous 
figures with the yellow-shaded histogram). Here we had to pick some variable which 
would characterize the true s 1 / 2 distribution. Two alternative choices which we also 
considered were the peak or the mean of the true s 1 / 2 distribution. All three of these 
variables are numerically quite close, with the peak value typically being the lowest, and 
the average value being the largest. In the end we chose (s 1 ^ 2 ) for its computational 
simplicity. This choice is rather inconsequential for our conclusions below, since we 
are introducing the {s 1 ^ 2 ) ave variable only for illustration purposes in Fig. ^. As we 
shall see, {s 1 ^ 2 ) ave actually cancels out in the final comparison between the next two 
variables. 



(s 1 / 2 )^ : this is the threshold of the true s 1//2 distribution, i.e. the minimum allowed 
value of s 1 / 2 . Since the minimum s 1//2 is obtained when the parent particles are produced 



-17 - 



0.2 



0.0 



a 0.5 



-1.0 



-1.5 1 — 1 
100 



gg-4j+;0 T 



(a) 



,1/2-1 



- As « (s ' (s , thr 

■Q 

O 

As'/ 2 =(s^(2m x )) peal -(s 1 / 2 ) lb 



As 1 / 2 = ( S ^ n (2m x )) peak -(s 1 /^) 



1/2n 



300 

m„ (GeV) 



1/2,1/2-, 



(c) 



.. As- = ( S ^) ave -(s^) thr 

O- .... _ 

As 1/2 = (sLf n (2m x )) peak -(s 1/2 ) tllr 



As 1 /^(sLf n (2m v )) peak -(s 1 / 2 ) a 



200 



300 

m x (GeV) 



400 



500 



0.2 



0.0 



S 0.5 



-1.0 



-1.5 1 — 1 
100 



gg-8j+£ T 



(b) 



, As 1 - (s'-)„ v( . (s"% r 

o- 

As 1 / 2 = ( S ^(2m x )) peak -( S '^) thr 

O Q. 

"'As 1 / 2 =(si/ 2 (2m x )) peak -(s 1 / 2 ) ave 



300 

(GeV) 



Xg-4j+£ T 
,1/2 / c l/2-i 



(d) 

,1/2-1 



.... As-=(s^) ave - (s V Z)tiir 
As 1 / 2 = (sL /2 (2m x )) peak -(s lG / 2 )' thr 



As 1 / 2 = (s^ n (2mJ) peak -( S 1 / 2 ) 



200 



300 

m x (GeV) 



400 



500 



Figure 8: Validity of the approximation ( |4.6p as a function of the LSP mass m x . The SUSY mass 
spectrum is fixed as in $3.16 ). In panels (a) and (h) we consider gluino pair production events, while 
in panels (c) and (d) we study associated gluino-LSP production. In panels (a) and (c) we force the 
gluino to decay to 2 jets as in (10.14 ), while in panels (b) and (d) each gluino decays to 4 jets as in ( 3.11 ). 
In each panel we compare the following three quantities: (s 1 ' 2 ) , which is the average of the true 

s 1 / 2 distribution; (s 1 / 2 ),, , which is the threshold of the true s 1 / 2 distribution; and ( s 1 J l ^ n (2m x ) ) , 

r V / peak 

1/2 

which is the location of the peak of the s' in (2m x ) distribution. 



at rest, (s 1 ' 2 )^ is nothing but the sum of the parent particle masses, as indicated in 
eqs. ( |4.1| ) and Q4.2|). Therefore, (s 1 ^ 2 ) i?tr is precisely the parameter that we would like 
to measure, in order to determine the true mass scale of the parent particles. 



^mm(^ m x)) • this is the parameter defined in eq. fl4.5|) , namely the location of 

/ peak 

1/2 

the peak of the •s n (j ri (2m x ) distribution, where we use the correct value for the invisible 
mass, in this case Mi nv = 2m x , since each SUSY event has two escaping LSPs. 



According to our empirically derived conjecture ( f4.6| ), the last two variables are approxi- 
mately equal, and the purpose of Fig. pi is to test this hypothesis, using the previously consid- 
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ered SUSY examples: gluino pair production (panels (a) and (b)), and associated gluino-LSP 
production (panels (c) and (d)). In panels (a) and (c) we force the gluino to decay to 2 jets 
as in ( 3.14 ), while in panels (b) and (d) each gluino decays to 4 jets as in ( 3.1 5| ) . Each line 
in Fig. E gives the fractional difference between a pair of |V2 

quantities as defined above. 

For normalisation we used the value of (s 1 ' 2 )^ , which is given by ( |4.1| ) for panels (a) and 
(b) and by (4.2) for panels (c) and (d). We vary the relevant part of the SUSY spectrum by 
changing the input value of the LSP mass m x and adjusting the other masses in accord with 

The main result in Fig. || is the comparison between the experimentally observable quan- 
tity (sllf n (2m x )] and the theoretical parameter (s 1//2 ) t/ir - As indicated by the red lines 

in Fig. H for the examples shown, those two quantities differ by no more than 10%, thus 
validating our conjecture (|4.6|) at the 10% level as well. We find this result quite intriguing. 
After all, we have not attempted any event reconstruction or decay chain identification, we 
are looking at very complex and challenging multijet signatures, and we have even included 
detector resolution effects. After all those detrimental factors, the possibility of making any 
kind of statement regarding the mass scale of the new physics at the level of 10% should be 
considered as rather impressive. 

We find it instructive to understand how we ended up with the observed precision, by 
comparing these two quantities |s^ n (2m x )j and (s 1 ^ 2 ) thr to the true s 1 / 2 as represented 

by its average (•5 1//2 ) troe - The blue lines in Fig. |8| show the fractional difference between 
{s 1/2 ) ave and (s 1/2 ) thr . We see that this difference varies by quite a lot, on the order of 
10-30% for gluino pair-production, but may get in excess of 150% for associated gluino-LSP 
production. As expected, (•s 1//2 ) croe is always larger than the threshold value (•s 1 '' 2 ),^, since 
the parent particles are typically produced with some boost, and the blue lines in Fig. |8| 
simply quantify the effect of this boost. 

On the other hand, the green lines in Fig. ^ represent the fractional difference (again 
normalised to (s 1 ^ 2 )^) between the measurable quantity {s l J^ n {2m x ) S j introduced earlier 

in eq. (|4-5| ), and the true energy scale of the events as given by (s 1 ^ 2 ) ave - We see that this 
time the fractional difference is negative, which simply reflects the fact that our variable 
^mm> being defined through a minimization condition, will always underestimate the true 
energy scale. The interesting fact is that while the blue and green curves in Fig. || have 
opposite signs, in absolute value they are very similar, leading to a fortuitous cancellation. 
The resulting discrepancy indicated by the red lines is therefore much smaller than either of 
the two individual errors indicated by the blue and green lines. 

It is now easy to understand qualitatively the origin of the approximate relation ( |4.6| ). 
Due to the boost at production, the true energy scale s 1 / 2 is larger than the threshold 
energy (•s 1//2 ) 4/ir by a certain amount. Later on, when we approximate s 1 / 2 with s^J 2 ^ we 
underestimate the true energy scale s 1 / 2 by more or less the same amount, bringing us 
back near the threshold (s 1//2 ) t/ir ,- As a result, the s l J 2 n distribution peaks very near the mass 
threshold (•s 1 ^ 2 ) t/tr which we are trying to measure in the first place. Of course, the proximity 
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of the s 1 J^ n peak to the threshold (•s 1//2 ) t?ir will be process dependent, but according to the 
examples considered here, holds to a remarkable accuracy. 

-1/2 

5. Correlation of the s^ in peak with the heavy particle mass threshold 



In the absence of a rigorous mathematical derivation, eq. ( [4.6[ ) should be considered simply 
as a conjecture. Nevertheless, once eq. Q4.6| ) is assumed to be approximately true, it allows 
us to measure the mass scale of the parent particles in terms of the hypothesized test mass 
rh x of the lightest invisible particle, e.g. the LSP in SUSY. For example, in the case of gluino 



pair-production in SUSY, we can use eqs. (4JJ and (4.6) to obtain a measurement of the 
gluino mass 

rhg(m x ) « - (s^ n (2m x )J (5.1) 

Z \ / peak 

as a function of the trial LSP mass rh x . Similarly, we can measure the gluino mass even in 
the much more challenging case of associated gluino-LSP production: from eqs. (|4.2| ) and 
( |4.6| ), we obtain 

rhg(m x ) « {s]£ n {2m x )) (5.2) 

V / peak 

As evidenced from eqs. (O) and (O), these measurements are very straightforward, since 



the only experimental input needed for them is the location of the peak of our all-inclusive 

1/2 

global variable s^ in . One should not be bothered by the fact that we did not get an absolute 
measurement of the gluino mass, but only obtain it as a function of the LSP mass. This is a 
well-known drawback of the other common mass measurement methods as well. For example, 
the classic M?2 endpoint analysis only yields the heavier parent mass as a function of the 
lighter child mass [9]. Similarly, the measurement of a single endpoint in some observable 
invariant mass distribution provides only a single functional relation between the masses of 
the intermediate particles in the decay chain, and by itself does not measure the absolute 
scale. In this sense, our measurement ( |5.lD is on equal footing with the more traditional 
methods. 

However, it is worth emphasizing the advantage of our method in the case of asymmetric 
events, where the parent particles are very different. An extreme version of such events is 
provided by the associated gluino-LSP production considered earlier. Under those circum- 
stances, the standard Mt2 method does not apply, while the single decay chain in the event 
may prove to be too short or too messy to provide a clean measurement through the invariant 

-1/2 

mass endpoint method. In contrast, we can still utilize s '■ for the measurement indicated 



in (|5.2| ) and a corresponding gluino mass determination. 

Let us now see how well the proposed measurements (^^) and (|5.2| ) will do for each of 
the SUSY examples considered in the previous section. In Fig. ||(a) we used eq. ( |5.1|) to 
convert our previous measurements of the various s^ n (2%) peaks in Figs. || and ^ into a 
corresponding gluino mass measurement. The red (blue) dashed lines correspond to the case 
of 4-jet (2-jet) gluino decays as in (|3l5D ( (EI3D ). We show results for the same four study 
points used in the four panels of Figs. || and [f| and the open circles mark the locations of the 
true masses (m x ,rrig), for each study point. 
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Figure 9: The correlation between the test LSP mass rh x and the corresponding gluino mass rhg, 



derived from (a) our proposed measurement ( |5.l| ) in gluino pair-production events, or (b) our proposed 
measurement ( fij| ) in associated gluino-LSP production events. Red (blue) lines correspond to the 
case of gluino decays to 4 jets as in ( 3.1^) (gluino decays to 2 jets as in (3.14)). The black dotted lines 
in panel (a) indicate the theoretically derived correlation from an ideal Mti endpoint analysis, i.e. 
assuming perfect resolution of the jet combinatorial ambiguity and ignoring any detector smearing. 
The open circles mark the locations of the true masses (m x , mg), for each of our four study points. 



The quality of the measurement ( |5.1| ) can be judged from the proximity of the experi- 
mentally derived m g (m x ) curves shown in the figure to the exact location of the true masses 
(m x ,m y ). We see that both the red and blue curves in Fig. |9|(a) pass very close to the true 
answer, especially for the study points with lower m x . In fact, we obtain a better measure- 
ment from the more complex 8-jet events (the red curves). At first sight, this may seem 
counterintuitive, until one realizes that the more visible objects are present in the event, the 
smaller the effect of the missing particles, and hence the smaller the error due to our approx- 



imation ( |2.4| , 2.5). Such multijet events appear very challenging to be tackled by any other 

means. For the sake of comparison, the black dotted lines in Fig. |9|(a) show the theoretically 

derived correlation from an ideal Mt2 endpoint analysis, i.e. assuming perfect resolution of 

the jet combinatorial ambiguity and ignoring any detector resolution effects. Comparing the 

red line from our measurement ( pM| ) to the ideal Mt2 line, we are tempted to conclude that, 
„i/2 — 

in essence, our s^ iri variable contains pretty much the same amount of information as Mt2- 
The big advantage of s^ iri , however, is the fact that we can obtain this information at a much 
lower cost in terms of analysis effort. 

Finally, in Fig. ||(b) we show our results from the analogous measurement ( ^.2[ ) in the 
case of associated gluino-LSP production. Here we also consider two different options for the 
gluino decay — 2 jet decays as in ( |3.14 ) (blue lines), or 4 jet decays as in ( ^T5|) (red lines). 
We then plot the resulting functional dependence rhg(rh x ) for each of the four study points 
considered earlier. Comparing Fig. |9|(b) to Fig. ^(a) which we just discussed, we arrive at 
very similar conclusions: the measurement ( |5.2[ ) is still quite accurate, and the superior result 
is provided by the more complex topology. Notice that here we do not show any Mr2-based 
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results, since the concept of Mt2 can not be applied to an extremely asymmetric topology 
like this one. 

6. The impact of initial state radiation and multiple parton interactions 

Up to now we have been discussing the ys variable of the primary parton-level hard scat- 
tering (HS). In principle, vs can be measured exactly, whenever we could both detect and 
identify the decay products of the heavy particles which were initially produced in the HS. 
Unfortunately, in reality it is rather difficult to measure V§ directly, for a couple of reasons: 

1. Omitting relevant particles from the y§ calculation. This case arises whenever some 
of the decay products resulting from the HS are not detected. For example, this may 
happen due to the imperfect hermeticity of the detector, where some of the relevant 
decay products are lost down the beam pipe. Fortunately, in reality this effect is pretty 
small. A much more serious problem arises whenever there are invisible particles Xi 
(see Fig. [l]) among the relevant decay products. Then, a relatively large fraction of the 
initial \f§ may go undetected, as can be seen by comparing the V§ m i n distributions in 
Figs. ||f|7| to the respective true (yellow-shaded) V§ distributions. 

2. Including irrelevant particles in the \f§ calculation. In general, any given event will 
contain a certain number of particles which will be seen in the detector, but did not 
originate from the primary HS. Initial state radiation (ISR), multiple parton interactions 
(MPI) and pile-up are the main examples of processes contributing to this effect. The 
pile- up effect can be controlled by a suitable Az cut, removing from consideration tracks 
which do not appear to originate from the primary vertex. However, ISR and MPI can 
be a serious problem. Including the extra particles will necessarily lead to an increase 
in the measured value of y§. In order to emphasize this difference, in the rest of this 
section we shall be using a prime to designate the experimentally measured quantities 
which include the full ISR and MPI effects (vF and v¥ m i n , correspondingly). 

Our proposal for dealing with the first of these two problems was to introduce the V§ m in 
variable in lieu of the true \f§. We then found an interesting empirical correlation (|4.6| ) 
between yfi m i n {2m x ) and the new physics mass scale. Now we shall turn our attention to 
dealing with the second problem, namely the fact that m in > V$min- 

Before we begin, we should mention that, depending on the particular circumstances 
and/or the goal of the experimenter, there may be certain situations where the inequality 
Vs^min > \f§min may not represent an actual problem. For example, if one is simply trying 
to measure the total energy in the observed events and not just the energy of the HS, then for 
missing energy events the relevant quantity of interest would be V^ m j„ itself, which would 
still be given by the expression ( |1.4| ) derived in Section |2|. There may also be situations 
where the ISR and/or MPI products may be reliably identified and excluded from the V§ m i n 
calculation. For example, consider a lepton collider and a missing energy signature with any 
number of jets and/or leptons. Since MPI is absent, while ISR and beamstrahlung would 
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only contribute photons, there will be no confusion with regards to which particles are due 
to ISR and which are coming from the HS. The analogous example at hadron colliders would 
be a signature containing anything but QCD jets. In what follows we shall ignore such trivial 
cases and instead focus on the much more challenging case of hadron colliders and jetty 
signatures, where the ISR/MPI products cannot be easily recognized. 

In the absence of any reliable methods for resolving the jet combinatorial problem on an 
event by event basis, one is left with two options. First, one may try to compensate for the 
ISR/MPI effects on the global VS'min distribution. In order to do this, one needs to know 
how ISR/MPI would affect the original V§ m in distribution. Ideally, this information should 
be measured from real data, using some Standard Model process as a standard candle. For 
example, Drell-Yan can provide the relevant information for a qq initial state [59], while tt 
can be used to study the gg initial state. Alternatively, one may calculate the ISR effects 
from first principles in QCD. Both of these approaches will be pursued in a future work [60]. 

A second approach would be to design and apply cuts which would minimize the ISR 
and MPI effects on the calculation of V^mm- Unfortunately, this is rather difficult to do in 
a model-independent fashion, since the size of the ISR effect is very model-dependent and 
depends on many factors: the energy of the collider (Tevatron or LHC), the mass of the 
produced particles, the identity of the partons initiating the HS, etc. Therefore, the optimal 
method to compensate for the ISR effect will also depend on all of these factors and will need 
to be decided on a case by case basis. 

For the purposes of the current study, we shall use a simple cut-based approach as 
discussed here, postponing the more complete treatment for [60]. To this end, we need to 
identify some global property of the ISR and MPI products which would distinguish them 
from the HS. Since it is well known that ISR and MPI peak in the forward region, it is 
natural to consider the pseudorapidity 77 as a simple cut variable. The energy distributions 
as a function of \rj\, for a few representative cases are shown in Fig. 10. We again consider 
the processes of gluino pair production (Figs. [U](a) and 0(b)) and associated gluino-LSP 
production (Figs. [0(c) and 0(d)). In each case, the gluino decays to 2 jets as in ( 3. 14Q . We 
choose to show the two extreme cases for the mass spectrum considered earlier: m x = 100 
GeV (Figs. |T|(a) and 0(c)) and m x = 400 GeV (Figs. 0(b) and 0(d)). The gluino mass is 
still fixed according to the gaugino unification relation ( |3.16| ). The black histograms in Fig. 10 



represent our previous results from Section |4| without any ISR or MPI effects, while the green 
(red) histograms include the effect of ISR (MPI). Finally, the blue histograms include both 
the ISR and MPI effects. The plots in Fig. [0 are normalized as follows. For each event, say 
the i-th one, we add the energy deposits in all calorimeter towers at a given \r)\, then divide 
the sum by the total energy Ei observed in the i-th event and the total number of events N, 
and finally enter the result into the corresponding \rj\ bin. It is easy to see that this ensures 
that the final distributions are unit-normalized. 



Fig. 10 shows that, as expected, the ISR and MPI effects appear mostly in the forward 
region. Therefore, by applying a simple |^| < r\ max cut, we could reduce their impact. Of 
course any such rapidity cut would essentially bring us back closer to the transverse quantities 
from which we were trying to escape from the very beginning. Furthermore, such a simple- 
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Figure 10: Energy distributions as a function of \rj\, for some of the SUSY examples considered 
earlier: 2-jet gluino decays from gluino pair production with (a) m x = 100 GeV or (b) m x = 400 
GeV; and from associated gluino-LSP production with (c) m x = 100 GeV or (d) m x — 400 GeV. The 
color scheme is such that the black histograms correspond to our previous results from Section ^ in 
the idealised case without ISR or MPI, while the green (red, blue) histograms include the effect of 
ISR (MPI, both ISR and MPI). Here Ei is the total energy measured in the i-th event, and N is the 
total number of events. As a result, all distributions shown in the figure are unit-normalized. 



minded procedure would introduce an uncontrollable systematic error, which would have to 
be estimated on a case by case basis. For example, Fig. |i~0|(b) shows that when the spectrum 
is rather heavy, the ISR/MPI effects are relatively small and can probably be safely neglected 
altogether, while Figs. |lO|(a), [H](c) and 0(d) reveal a significant ISR/MPI pollution for a light 
SUSY spectrum. One should also keep in mind that our conjecture (|4.6|) is already subject 
to a certain systematic error, whose size sets the benchmark for the ISR/MPI elimination 
study. With those caveats, we choose our cut at r] max = 1.4, which is nothing but the end of 
the barrel and beginning of HE/HF calorimeters in CMS. This choice makes good sense from 
an experimentalist's point of view, since the segmentation and performance of the HE/HF 
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Figure 11: Unit-normalized distributions of y§min{^ m x) an< ^ v / ^'mm(2m x ) for the SUSY examples 
considered in Fig. [!(]. The color scheme is the same as in Fig. |T(]. The blue histograms include both 
the ISR and MPI effects, and represent the actually measured v / i' rn i„(2m x ), while the green (red) 
histograms include the effects of ISR (MPI) only. All three of those distributions are subject to the 
\r)\ < 1.4 cut discussed in the text. For comparison, we also show our previous results from Section 
corresponding to the HS only (without any ISR or MPI effects) and without an n cut. In particular, 
the black solid histograms in Fig. [TI] represent our previous results for the quantity y/I m i n (2m x ), 
while the black dotted (yellow-shaded) histograms give the true s 1 / 2 distribution, whose threshold is 
the parameter to be measured. 



calorimeters are relatively worse to begin with. 

Let us now revisit some of the V»min distributions from Section || and incorporate suc- 
cessively the effects of ISR and/or MPI. Fig. 11 shows our results for the same four SUSY 
examples from Fig. 1C. The green (red) histograms include the effect of ISR (MPI) alone, 
while the blue histograms include both the ISR and MPI effects, and thus represent the true 
measured quantity y / i' m j n (2m x ). All three of those distributions are subject to our |^| < 1.4 
cut. For comparison, we also show our previous results from Section |], corresponding to the 
HS only (without any ISR or MPI effects) and without an 77 cut. In particular, the black 
solid histograms in Fig. 11 represent our previous results for the Quantity v<5min 

(2m x ), while 
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the black dotted (yellow-shaded) histograms give the true s 1//2 distribution, whose threshold 
is the parameter that ideally we would like to measure. 

Fig. 11 confirms that the ISR and MPI effects shift the original HS distribution \Th m in 



(black histograms) into a harder ys' m j n distribution (blue histograms), even after applying 
the rj cut. The size of this effect depends on the mass spectrum: it is more pronounced 
when the spectrum is light 5 , as in Figs. |TT|(a), [Tl|(c) and [TT|(d) . In the worst case scenario of 
Fig. |ll](c) the location of the s l J^ n peak shifts by almost a factor of two. On the other hand, 
for the best-case scenario of Fig. 11(b) the shift is rather small. By comparing the green and 
red histograms, we can also deduce the relative importance of ISR versus MPI. We see that 
the two effects are roughly comparable in size, but as a rule, the red histograms are shifted 
further along, which suggests that MPI has a somewhat higher impact than ISR, indicating 
the importance of understanding the full structure of the underlying event at the LHC. The 
general conclusion from Fig. 11 is that our mass measurement method proposed in Section ||] 
is likely to work much better if the new particle spectrum happens to be relatively heavy. 
This assumption is not unreasonable: if the new physics spectrum were too light, then it 
might have already been ruled out directly or indirectly, and if not, then due to the higher 
production cross-sections, there should at least be sufficient statistics to attempt some sort of 
exclusive reconstruction. In this sense, for the case where \/~§ m in is most likely to be useful, 
ISR and MPI are least likely to be a problem. 

We are now in a position to repeat our mass measurement analysis from Section || with 
the inclusion of ISR and MPI, while ignoring the forward calorimetry through an \rj\ < 1.4 
cut. Our results are shown in Fig. 12. Comparing Figs. || and 12, we see that the inclusion 
of ISR/MPI deteriorates the mass measurement, most notably for light SUSY mass spectra 
with m x ~ 100 — 200 GeV. This should not be surprising, given what we have already seen in 
Figs. [l0| and 11. Nevertheless, for heavier SUSY spectra the precision remains relatively good, 
typically on the order of 10%, even for the most challenging cases of associated gluino-LSP 
production. 



7. Summary and conclusions 

Anticipating that an early (late) discovery of a missing energy signal at the LHC (Tevatron) 
may involve a signal topology which is too complex for a successful and immediate exclusive 

~l/2 

event reconstruction, we proposed a new global and inclusive variable s'-, defined as follows: 
it is the minimum required center-of-mass energy, given the measured values of the total 
calorimeter energy E, total visible momentum P, and/or missing transverse energy $t in the 
event. Our variable has several desirable features: 

• It is global in the sense that it uses all of the available information in the event and not 
just transverse quantities, for example. 

J Notice that for a given value of m x , the relevant mass scale 2m g = 12m x in Figs. |il|(a ) and 0(b) is 
almost twice as large as the corresponding mass scale trig + m-o = 7m x in Figs. |ll|(c) and |ll|(d). 
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Figure 12: The same as Fig. but including the effects of ISR (green), MPI (red), both ISR and 
MPI (blue). The variable s'^(2m x ) here is 
coverage of the CMS barrel calorimeter only. 



1/2 

MPI (blue). The variable s' min (2m x ) here is calculated with a cut of \n\ < 1.4, corresponding to the 



• It is inclusive in the sense that it does not depend on the specific production process, 
or particular decay chain. Consequently, it is also very model-independent and does 
not require any exclusive event reconstruction, which may be a great advantage in the 
early days of the LHC. 

• It is theoretically well defined and as such has a clear physical meaning: it gives the 
minimum total energy which is consistent with a given observed event. This intuitively 
clear physical picture allowed us to correlate it with the mass threshold of the new 
particles as in eq. (f4.6|), which turned out to work surprisingly well. In contrast, it is 
generally difficult to correlate a bump in a purely transverse quantity like Et or Ht to 
any physical mass parameter in a model-independent fashion. 

In Section ^ we derived a simple formula for s^ in in terms of the measured E, P z 
and Et- The formula is in fact completely general, and is valid for any generic event shown 
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in Fig. |l[ with an arbitrary number and/or types of missing particles. Therefore, it can be 
applied equally successfully to SM as well as BSM missing energy signals. 

In Sections || and |] we identified two useful properties of the s 1 J^ n variable. First, its shape 
matches the true s 1 / 2 distribution better than any of the other global inclusive quantities 
which are commonly discussed in the literature. More importantly, when we create the 
s 1 ^ 2 (Mi nv ) distribution with the true value of the invisible mass Mi nv , its peak is very close 
to the mass threshold of the parent particles originally produced in the event. This conjecture, 
summarized in eq. (^1]) , allows us to obtain a rough estimate of the new physics mass scale, as 
a function of the single parameter Mj n „. For example, in i?-parity conserving super symmetry, 
where Mi nv = 2m x , we derive a relation between the heavy superpartner mass and the mass 
of the LSP, as shown in Fig. ^. 

1/2 

Before we conclude, we should comment on several other potential uses of the s^ in vari- 

1/2 

able. Before we even get to the discovery stage, s^ n (0) can already be used for background 
rejection and increasing signal to noise, just like Mt2(0) [25]. In particular, it is interesting 

1/2 

to explore the correlations between s^ in and the other global inclusive variables discussed in 
Section | [61]. While we did not include any SM backgrounds in our SUSY plots, we expect 
that the presence of SM backgrounds will not affect either the existence or the location of 

1/2 1/2 

the new physics 8^(0) peak. At large values of ^^(0), where a new physics signal is most 
likely to appear, any SM background will be rather smooth and featureless, so that it can be 
safely subtracted away through a side-band method. 

Another possible application of sllf n {0) is at the trigger level. In Section || we already 

1/2 

saw that s^ n (0) is superior to both Ht and $t in identifying the scale of the hard scattering. 
At the same time, there exist dedicated Ht and $t triggers, motivated by the sensitivity of 
those variables to the relevant energy scale. Given that our variable is doing an even better 

1/2 

job in this respect, we believe that the implementation of a high-level s^ in (0) trigger should 
be given a serious consideration. 

-1/2 

As we have been emphasizing throughout, a major advantage of s'- is that it does not 
require any explicit event reconstruction and thus it is very model- independent. We should 
mention that to some extent, these properties are also shared by the MTGen variable proposed 
in [33]. In calculating MTGen-, one considers all possible partitions of the visible particles Xi 
in the event, thus effectively eliminating the model-dependence which stems from assuming a 

1/2 1/2 

particular topology. While MTGen and s^ in are similar in this respect, we believe that s^ in 

1/2 

has three definite advantages — first, it is much, much easier to construct. Second, s^ in 
can be applied to extreme asymmetric topologies where the second side of the event yields 
no visible particles. A simple example of this sort was the associated gluino-LSP production 
considered in Figs. [|, [7], |8|(c,d), ^(b) and p^(c,d). Finally, the interpretation of sHf n involves 
reading off a peak, while MTGen requires reading off an endpoint. The former is much easier 
than the latter: for example, a peak would still be recognizable in the presence of large 
backgrounds. In contrast, an MTGen endpoint can fade out due to a number of reasons, 
including detector resolution, combinatorial background, etc. On the other hand, MTGen 
(and more generally, the Mt2 class of variables) is better behaved in the presence of ISR. 
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More specifically, the endpoints of the MxGen and Mtz distributions in general do shift in 
the presence of ISR, and their explicit dependence on the "upstream" transverse momentum 
has to be calculated on a case by case basis [54]. However, the nice feature of both MTGen 
and Mt2 is that when the test mass m x becomes equal to its true value m x , there is no such 
shift and the endpoint remains intact even in the presence of arbitrary ISR. In contrast, as 
discussed in Section ^, s 1 J^ n is always affected by ISR to some extent, requiring some sort of 
correction. 

-1/2 

In conclusion, we reiterate that perhaps the most important advantage of s^ in is that it 

-1/2 

is readily available from day one. We are therefore eagerly looking forward to the first s' in 
plots produced with real LHC data. 
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